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ABSTRACT. We study small random perturbations by additive white-noise of a spatial discretization 
of a reaction-diffusion equation with a stable equilibrium and solutions that blow up in finite time. We 
prove that the perturbed system blows up with total probability and establish its order of magnitude 
and asymptotic distribution. For initial data in the domain of explosion we prove that the explosion 
time converges to the deterministic one while for initial data in the domain of attraction of the stable 
equilibrium we show that the system exhibits metastable behavior. 



1. Introduction 
We consider small random perturbations of the following ODE 

r u> = 2(_tf 1+ tf 2 ), 

(1.1) < VI = UUi+i-lUi + Ui-x) 2<i<d-\, 

[ U' d = %(-U d + U d -i+hg(U d )). 

Here g : R — » R is a reaction term given by g(x) = (x + ) p — x with p > 1, and h > is a parameter. 
We also impose an initial condition Uq € M. d . This kind of systems arise as spatial discretizations of 
diffusion equations with nonlinear boundary conditions of Neumann type. In fact, it is known that as 
h — y solutions to this system converge to solutions of the PDE 

= u xx (t,x) 0<x<l,0<t <T, 
= < t < T, 

= 0<t<T, 

= Uq{x) < X < 1. 

This and more general reaction-diffusion problems including for instance the possibility of a non- 
linear source term like g and other type of boundary conditions appear in several branches of pure 
and applied mathematics. They have been used to model heat transfer, exothermic chemical reactions, 
population growth models, geometric flows, etc. 

An important feature of this type of problems is that they admit solutions which are local in time, 
with the possibility of blow-up in finite time. The asymptotic behavior of solutions to (11.11 ) can 
be briefly summarized as follows (we give a detailed description afterwards): the system has two 
equilibriums Uq = and Uq = 1 . The first one is stable while the second is unstable. Hence, there 
exists a domain of attraction Do for the zero solution such that if Uq € Do then the solution U (?) = 
(U\(t),. . . ,U d (t)) with initial condition Uq is globally defined and U(t) —> as t — s> °°. There exists 
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also a stable manifold for the unstable equilibrium which is of co-dimension one and coincides with 
the boundary of Do- For Uq G D c the solution U blows up in finite time T = T(Uq). 

Since mathematical models are not exact, it is important to understand what changes arise in the 
behavior of the system when it is subject to perturbations. We study random perturbations given by 
additive white-noise. More precisely, we consider Stochastic Differential Equations (SDE) of the form 

( dU l ; £ = l(-U^ + U^)dt + edW u 

(1.2) < dU- l ' e = jr(U"_l\ -2U"' e + U.'j_])dt + edWi 2<i<d-l, 
{ dU^ £ = ^(-U^ + U^+hgiU^dt + edWj, 

which can be written in matrix form as 

(1.3) dU luE = (-A£/"' E + p(U%' e )e d )dt + edW. 

Here W = (W\ Wd) is a d— dimensional standard Brownian motion, £ > is a small parameter and 
e d = (0, . . . , 1) is the d-th canonical vector on W 1 . In the sequel we use U u,£ for a solution to (11.21) 
with initial condition U" ,E (0) = u 6l . In the case £ = we are left with the deterministic equation 
and so we use the notation U" := U u '° to denote a solution to (ll.ll ). 

The field b(U) := -AU + \g(U d ) e d is a gradient (b = — V(j)) with potential given by 

1,. 2/ \U+\P +l Uf- 



ty(U) = -(AU,U)-- 



2 X ' ' h\ p + l 2 

The SDE associated to this energy functional can be compared with the classic double-well potential 
model, which we now briefly summarize. We refer to iflOl p. 294] for a more detailed description. 

In the double-well potential model one considers a stochastic differential equation of the form 

(1.4) dX e = r(X e )dt + edW 

where W is a standard <i-dimensional Brownian motion and r is a globally Lipschitz gradient field 
over M d given by the double-well potential (j). More precisely, this potential (j) possesses exactly 
three critical points: two local minima p and q of different depth and a saddle point z with higher 
energy, that is (j)(z) > §(p) > §(q) . Each minimum corresponds to a stable equilibrium and hence 
for initial data lying outside the stable manifold of z, the deterministic system (£ = 0) converges to 
one of them depending on the initial condition. When considering random perturbations, for compact 
time intervals the stochastic system converges as £ — > to the deterministic one uniformly but the 
qualitative behavior of the perturbed system is quite different from that of the deterministic solution 
for large times. If the potential grows fast enough at infinity the resulting stochastic system admits 
a stationary probability measure which converges to a Dirac delta concentrated at the bottom of the 
deepest well q. Hence, for initial data in the domain of attraction of the shortest well p we observe 
that 

(i) Due to the action of the field r, the process is attracted towards the bottom of the shortest well 
p; once near p, the field becomes negligible and the process is then pushed away from the 
bottom of the well by noise. Being apart from p, noise becomes oveipowered by the field r 
and this allows for the previous pattern to repeat itself: a large number of attempts to escape 
from the given well, followed by a strong attraction towards its bottom. This phase is known 
as thermalization. 

(ii) Eventually, after many frustrated attempts, the process succeeds in overcoming the barrier of 
potential and reaches the deepest well. Since the probability of such an event is small, we 
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expect this tunneling time to be exponentially large. Moreover, due to the large number of 
attempts that are necessary, we expect this time to show little memory, 
(iii) Once in the deepest well, the process behaves as in (i). Since the new barrier of potential is 
higher, the next tunneling time is expected to happen on a larger time scale. 

This description was proved rigorously in (3l |5j d 01 [H using different techniques. The phenomenon 
is known as metastability . For a detailed description of it we refer to |[T0l . 

Coming back to our potential (j), the situation is slightly more complex. Instead of having a deepest 
well, we have a direction along which the potential goes to — oo and, hence, the size of the "deepest 
well" is now infinity and there is no return from there. Moreover, since the potential behaves like 
— s p+l in this direction, if the system falls in this "well", it reaches infinity in finite time (explosion). 

The purpose of this paper is to study the metastability phenomenon for this kind of potentials 
where there is a shortest (finite) well and a deepest well which yields to infinity in finite time. The 
ideas developed here can be extended to other systems with the same structure. The typical situation 
with this kind of geometry is the case of reaction-diffusion equations where the reaction comes from 
a nonlinear source with superlinear behavior at infinity such as 

i p 

u t — Uxx ~r M + , 

with p > 1, in a bounded domain of M and homogeneous Dirichlet boundary conditions. In this case 
the diffusive term pushes the solution towards zero (a stable equilibrium) while the source u p + pushes 
it to infinity. In this situation we expect the same behavior as the one of solutions to (11.2I ). 

Since the drift in (11.21) is not globally Lipschitz, we are only able to prove the existence of local 
solutions and in fact, explosions occur for solutions of (1 1 -2b - In particular, classical large deviation 
principles as well as other Freidlin-Wentzell estimates do not apply directly. All of these results deal 
with globally Lipschitz coefficients. Also, the loss of memory for the tunneling time was proved 
only in the globally Lipschitz case where explosions do not occur. The only exception is the work 
of Azencott [2] where locally Lipschitz coefficients are considered and explosions are allowed, but 
the large deviations estimates developed there apply only to neighborhoods of solutions which do not 
explode in a fixed time interval (and hence the perturbed system is automatically defined in the whole 
interval for £ small enough). In that work the author also considers the exit from a domain problem, 
but explosions are not allowed in his analysis. 

As opposed to this last case, we specifically focus on trajectories that blow up in finite time. The 
asymptotic behavior (as £ — > 0) of the explosion time for (11.21) is not understood yet, and this is the 
goal of this article. 

In order to study this kind of systems, localization techniques may be applied but this has to be done 
carefully. The main difficulties lie in (i) the geometry of the potential (and its respective truncations) 
which is far from being as simple as in the double-well potential and (ii) the explosion phenomena 
itself. Localization techniques apply reasonably well to deal with the process until it escapes any 
bounded domain, but dealing with process from there up to the explosion time requires different tools, 
which include a careful study of the blow-up phenomenon. Clearly, localization arguments are useless 
for this last part. 

The paper is organized as follows. In Section [2] we give the necessary definitions, review some 
Freidlin-Wentzell estimates and detail the results of this article. Section [3] is devoted to giving a 
detailed description of the deterministic system (ll.ll ). In Section @] we begin our analysis of the 
stochastic system. We prove that explosions occur with probability one for every initial datum. In 
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Section [5] we prove that for initial data in the domain of explosion, the explosion time converges to 
the deterministic one as 8 — > 0. Finally, throughout Section [6] we study the characteristics associated 
to metastability for initial datum in the domain of attraction of the origin: exponential magnitude of 
the explosion time and asymptotic loss of memory. 



2. Definitions and results 

2.1. Solutions up to an explosion time. Throughout the paper we study stochastic differential equa- 
tions of the form 

(2.1) dX = b(X)dt + edW 

where £ > and b : R d — > W 1 is locally Lipschitz. It is possible that such equations do not admit 
strong solutions in the usual sense as these may not be globally defined but defined up to an explosion 
time instead. We now formalize the idea of explosion and properly define the concept of solutions for 
this kind of equations. We follow Q- 



Definition 2.1. A solution up to an explosion time of the stochastic differential equation (12.11 ) on the 
probability space (£2, J ,P), with respect to a filtration (j t )t>o satisfying the usual conditions and a 
fixed Brownian motion {W t: f t )t>o with (a.s. finite) initial condition £, is an adapted process X with 
continuous paths taking values in W 1 U {°°} which satisfies the following properties: 



If we define x" = inf{f > : \X(t) \ = n} then for every n > 1 we have 

ffAT" 

and 




\b(X(s))\ds<+°° =1 V0<?<+oo 



p(x{tAi n ) + £b(X(s))t {s < t , } ds + eW(t Ax"); V < t < +<~J =1. 

• X has the strong Markov property, i.e. if we note x := lim„^ +0 oX" and x is a stopping time 
of X then, conditional on x < X and X(x) = x, the future {X + (t) = X(t + x) : t < x - x} is 
independent of the past {X(s) : s < x} and identical in law to the process started at x. 

We call x the explosion time for X. Notice that the assumption of continuity of X in W 1 U {°°} 
implies that 

x = inf{? >0:X(t) (£R d } and X(t-) = X(%) = °o on {x < +°o}. 

We stipulate that X(t) = oo provided that x < t < +oo but we do not assume that \im t ^ +C oX(t) exists 
when x = +oo. 

Notice that the assumption of finiteness of ^ grants us P(x > 0) = 1. Also, if P(x = +oo) = 1 then 
we are left with the usual definition of strong solution to the equation. 



Remark 2.1. It can be proved that if b € C 1 (R ) then there exists a unique solution of (12.11 ) up to an 
explosion time (see (SHH). 
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2.2. Freidlin-Wentzell estimates. One of the most valuable tools in the study of perturbations by 
additive white noise of an ODE is the Freidlin-Wentzell theory, whose main results we briefly describe 
here. 

LetX*> 8 be a solution to the SDE 

dX x ' E = b(X x *)dt + edW 

with initial condition x € M. d , where b is globally Lipschitz with Lipschitz constant K. Fix T > and 
let Px' T denote the law of X x,e on C([0,T],M. d ). Let us also consider X x the unique solution to the 
deterministic equation 

X(t) = ~b(X(t)) 

with initial condition x € W 1 . 

Theorem 2.2 (Freidlin and Wentzell, (31). For each x € W 1 and T > the family (Px' T )e>o satisfies a 
large deviations principle on C([0, T],M. ) with scaling £~ 2 and (good) rate function Ij given by 

{\ So \§i s ) ~ ^( ( P( S )) | 2 ds if cp is absolutely continuous and (p(0) = x 
+oo otherwise 

As a matter of fact, we need only the following weaker statement for our analysis: for every fixed 
T > and 8 > there exist positive constants C\ and C2 depending on T, 8 and K such that for all 
0<£< 1 

(2.2) supPj sup \^ X ' Z {t)-X x {t)\ >5j <Cie"&. 

.reR rf \te[o,r] / 

2.3. Main results. We now state the main results of the article. The first of them concerns the explo- 
sion time of solutions to 111. 21 . In the following P u denotes the law of the solution to (11.21 ) up to the 
explosion time x" with initial condition u. When the initial condition is clear we often write x e instead 
of x" to simplify the notation. 

Theorem 2.3. Let U"' e be a solution to Oil . Then P i( (x E < 00) = 1. 

Let us notice that this result establishes a first difference in behavior with respect to the deterministic 
system. While global solutions exist in the deterministic equation, they do not for the stochastic one. 

We then focus on establishing the order of magnitude and asymptotic distribution of the explosion 
time for the different initial conditions u € R . We deal first with initial conditions in the domain of 
explosion D e and show the following result. 

Theorem 2.4. Given 8 > and u G D e we have 

(2.3) limP H (|T e -T | >8) = 0. 

e->0 

Moreover, the convergence is exponentially fast. 

This last theorem shows that for small £ > the behavior of the stochastic system does not differ 
significantly from the deterministic one for initial conditions in D e . However, this is not the case for 
initial data in the domain of attraction of the origin. Here is where important differences appear and 
where characteristics associated with metaestability are observed. In order to properly state the results 
achieved in this matter, we need to introduce some notation. 
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For each £ > we define 

p E = inf{f >0:P (x e >p e ) <e- 1 } 

which is well defined since Po(x e < +°°) — 1 f° r every £ > 0. We first show that the family (P E ) e >o 
verifies 

lim£ 2 logP e = A 

with A := 2((j)(l) — (j)(0)). In fact, we prove the stronger statement featured in the following theorem. 
Theorem 2.5. For each u 6 Do and 8 > 

/ A-S A+5 \ 

lim P u < x e < e~Z~ =1, 
where the convergence is uniform over compact subsets of Do. 



This theorem characterizes the asymptotic order of magnitude of the explosion time for any initial 
condition u € Dq. Regarding its distribution, we show the asymptotic loss of memory in our last result. 

Theorem 2.6. For each u G Do and t > 

limP l( (x E >?p E ) =e~' 

£H>0 

where the convergence is uniform over compact subsets of Do. 

3. The deterministic system 

Throughout this section we state some properties and study the behavior of solutions to (ll.lb . This is 
carried out in (TJ for solutions with nonnegative initial conditions. The purpose of this section is to 
extend the analysis in [TJ to any arbitrary initial data u € M. d . 

Let us start by noticing that equation (11.11 ) can be written as 

U{t)=b{U{t)) 

for b = — V(j) where (]) is defined as 

1 ? / \TJ + \P +l JjJ 2 



2 X ' ' h\ p + l 2 

Here A is as in (11.2I )- (11.3I ). Notice that the potential (j) has exactly two critical points: 1 := (1, . . . , 1) 
and the origin. Both of them are hyperbolic. The origin is the only local minimum of (j) while 1 
is a saddle point. Our goal is to decompose W l into distinct regions, each of them having different 
asymptotic characteristics under our system. To be able to accomplish such decomposition we need a 
few results concerning solutions to dl.lb . We begin with the following proposition. 

Proposition 3.1. Let U = (U\,. . . ,Ud) be a solution to (ll.ll ). Then the application t i-> §(U (t)) is 
monotone decreasing. 

Proof. Since A is symmetric and U = —AU + y i g(Ud)e c i, a direct calculation shows that 
d -^P^- = (U(t),AU(t)) - lg(U d (t))U d (t) = -\U(t)\ 2 < 0. 

□ 
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Next we show that solutions to dl.ll ) satisfy a Maximum Principle. 
Lemma 3.2 (Maximum Principle). Let U = (U\, . . . , Ud) be a solution to dl.ll ). Then U satisfies 

(3.2) max \Uk(t)\ < max{ max |£4(0)|, max U c i(s)} 

k=l,...,d k=l,...,d 0<s<t 

Proof. We prove first that 

(3.3) max \U k (t)\ < max{ max \U k (0)\, max \U d (s)\} 

k=l,...,d k=l,...,d 0<s<t 

and then we check that if d3.3b holds then 

max{ max |£4(0)|, max |£/ f /(s)|} = max{ max |£4(0)|, max Ud (s)} 

k=l,...,d 0<s<t k=l,...,d 0<s<t 

which allows us to conclude (13.2b . Let j be the node that maximizes maxo< iV < f |t//(s)|. Let us observe 
that if j = d then (13.3b is immediately verified. Hence, we can assume that 1 < j < d. Consider 
to = min{Y G [0,?] : maxo<i< f | Uj (s) \ = \Uj(t')\}, the first time in which the maximum is attained. 
Note that |£//(/o)| = T^^k=i,...,d( max o<s<t l^it (*)!)■ If to = then 



max |I4(0)| > \Uj(t Q )\ = max max \U k (s)\ > max \U k (t)\ 

k=\,....d k=l,...,d \0<s<t / k=l,...,d 

and we get (13.3b . If to > we must consider two cases: Uj(to) > and Uj(to) < 0. If Uj(to) > then 
by definition of to we get that Uj(t ) > Uj(s) for all < s < t. From this it follows that Uj(t ) > 0. On 
the other hand, the choice of j guarantees that Uj(to) > Uk(to) for all k = 1, . . . ,d. This implies that 

u j< f o) = n(( t O+i( / o)-^o)) + (^-i(A))-^(A)))) <0 if Kj<d 



h 2 

and 



U{ (to) = Jr (U 2 (to) - U X (t Q )) < if j = 1. 



In any case we conclude that Uj(to) = and, in particular, that Uj + \(to) = Uj(to). We conclude 
that \Uj + i(to) \ = max / t = i i ... i( /(maxo<. v <r which allows us to repeat the same argument, now for 

7 + 1 instead of j. Thus, an inductive procedure eventually yields that Ud(to) = UAto). From here 
we obtain (13.3b if Uj(to) > 0. The case E//(/o) < is analogous. To conclude (13.2b we notice that if 
ti = mm{t' e [0,t] : max <. v < r \Ud(s) \ = \U d (t')\} > then Ud(ti) > because, otherwise, from (fTTTT) 
and (13.3b we get that U' d (t\) > which contradicts the definition oft\. □ 



As a consequence of the Maximum Principle we have the following characterization of globally 
defined solutions to dl.lb . 

Lemma 3.3. Let U be a globally defined solution to dl.lb . Then U is bounded. 

Proof. Let us suppose that U is not bounded. Then by the Maximum Principle we obtain that 

maxo<.s<( \U d (s) \ -> +°o as t -> +oo. 

1. Given M > we define ? M := inf{? > : \U d (t)\ > M}. From this definition it follows that 
\Ud(t M )\ > M and that \U d (t M )\ = max < s < (M \U d (s)\. If M > max^i,...^ |t4(0)| then ? M > and 
by the Maximum Principle we have \Ud-i (?m)| < U d (tM)- This gives us the inequality 

£/^m) > luZ(t M )-(^ + l)Ud(t M ). 



8 



PABLO GROISMAN AND SANTIAGO SAGLIETTI 



2. From here it is easy to see that if M is large enough we have that '■ [%, +°°) — > K is monotone 
increasing. This implies that for t > tM we have Ua(t) = maxo< s < ; \Ud(s) \ > M and, as a consequence, 
that U' d (t) > \U%(t) - + \)U d (t). If M is taken large enough then £/ verifies U' d (t) > \U%{t) for 
t >tM and, therefore, cannot be globally defined. This is a contradiction which implies that U must 
be bounded. □ 

From the previous lemma and the fact that (11.11) admits the Lyapunov functional (13.11 ) we obtain 
the following corollary. 

Corollary 3.4. Let U be a solution to (11.11) . Then either U explodes in finite time or is globally defined 
and converges to a stationary solution as t — > +oo. 

With this result at our disposal we can obtain the following theorem, whose proof is in ifTTl. 
Theorem 3.5. 

(1) Equation (11.11) has exactly two equilibriums U = and U = 1. The first one is stable and the 
second one is unstable. 

(2) Let u be a nonnegative initial datum such that U" is globally defined and lim t ^ +00 U u {t) = 1. 
Then 

• < v < u => U v is globally defined and lim U v (t) = 0. 

• u < v U v explodes in finite time. 

(3) Consider X > and a nonnegative initial condition u. Then there exists A- c > such that 

(a) X < X c U x " is globally defined and lim t ^ +oa U x " (t) = 

(b) X c < X U Xu explodes in finite time 

(c) X = X c U x " is globally defined and lim t ^ +oa U x " (t) = 1. 

This results allow us to give a good description of the behavior of the deterministic system U for 
the different initial conditions u € M. d . Indeed, we have a decomposition 

R d =D Q UWiUD e 

where Do denotes the stable manifold of the origin, W( is the stable manifold of 1 := (1, . . . , 1) and 
D e is the domain of explosion, i.e., if u € D e then U" explodes in finite time. The sets Do and D e are 
open in W 1 . The origin is an asymptotically stable equilibrium of the system. is a manifold of 
codimension one. Also 1 admits an unstable manifold of dimension one which we shall note by 'W^. 
This unstable manifold is contained in R^_, has nonempty intersection with both Do and D e and joins 1 
with the origin. An illustration of this decomposition is given in Figured] for the 2-dimensional case. 

4. Explosions in the stochastic model 

In this section we focus on proving that solutions to (11.21 ) blow up in finite time with probability 
one for any initial condition u 6 M. d and every e > 0. The idea is to show that, conditioned on non- 
explosion, the system is guaranteed to enter a specific region of space in which we can prove that 
explosion occurs with total probability. From this we can conclude that non-explosion must happen 
with zero probability. We do this by comparison with an adequate Ornstein-Uhlenbeck process. 

Proof of Theorem [PI Let Y y ' e be the solution to 

dY y e = - (AY y > e + \y y d ,Ze d) dt + edW 
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u = o 



Figure 1. The phase diagram of equation dl.lb . 



with initial condition Y y ' e (0) = y. Notice that the drift term is linear, and given by a negative definite 
matrix. Hence, Y y ' e is in fact a J-dimensional Ornstein-Uhlenbeck process which admits an invariant 
distribution supported in M. d . We also have convergence to this equilibrium measure for any initial 
distribution and therefore the hitting time of Y y,e of any open set is finite almost surely. 

On the other hand, since the drift term of ((4} is smaller or equal than b we can apply the stochastic 
comparison principle to obtain that U"- e (t) > Y y,e holds a.s. as long as U U,E is finite, if u > y. From 
here, the result follows applying the following lemma and the strong Markov property. □ 

Lemma 4.1. Consider the set 

e M := {y G R d : y k > for all < k < d - 1 , y d > M}, 

then we have 

lim inf P v (x e < °°) = 1. 

Proof. Consider the auxiliary process Z y,e := U y ' e — eW. Notice that this process verifies the random 
differential equation 

dZ y * = b(Z y * + eW)dt, Z- V ' E (0) = y. 

Let us also observe that Z V E has the same explosion time as U y,e . For each k G N let us define the set 
Au := {sup 0<f<1 |W ( /(f)| < k}. On we have thatZ^' 8 verifies the inequality 

(4.1) ^- > -AZ>'' E -^£e,- + -((Zf -zk%-Zf -ek)e d . 

Observe that (14.11 ) can be written as 

d7 y ' 8 

a J— > eZ .v.e + q + {z f _ zk y + > QZ y,* + qi 
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where Q G R dxd and q G R d both depend on e,h and k, but not on M. This allows us to conclude the 
inequality \Z y ' e \ < (M + |<?|)exp(|<2|) for all <t < 1. In particular, for the last coordinate we get 

>-aiM + a 2 (Zj e )f if0<?<l 
ZJ € (0)>M 

for constants 0Ci , CLj which do not depend on M. It is a straightforward calculation to check that solu- 
tions to this one-dimensional inequality blow-up in a finite time that converges to zero asM-> +oo. 
Therefore, for each k G N there exists such that P(/U) < inf ve0 M P v (x E < °°) for all M > M^. Since 
\imi c ^ +00 P(Ai c ) = 1, this concludes the proof. □ 



5. Convergence of t£ for initial conditions in D e 



This section is devoted to prove that for initial data in the domain of explosion of the deterministic 
system, the explosion time is of order one and, moreover, as £ — > converges to the explosion time of 
the deterministic system. Observe that do to the lack of boundedness this result do not follow from 
standard perturbation arguments for dynamical system (deterministic or stochastic). We first introduce 
the truncations of the drift that we use here to prove one of the bounds and we are going to make more 
profit of them in Section [6] when we deal with initial data in the domain of attraction of the origin. 



5. 1. Truncations of the potential and localization. The large deviations principle originally formu- 
lated by Freidlin and Wentzell for solutions of stochastic differential equations like (12.11 ) require a 
global Lipschitz condition on the drift term b. While this condition is met on the classic double-well 
potential model, it is not in our case. As a consequence, we cannot apply such estimates to our system 
directly. Nonetheless, the use of localization techniques helps us to solve this problem and allows us 
to take advantage of the theory developed by Freidlin and Wentzell despite the fact that our drift term 
is not globally Lipschitz. In the following lines we give details about the localization procedure to be 
employed in the study of our system. 

For every n G N let G n : R — > R be of class C 2 such that 

G n {u) = { Vrr-T rf|»l<» 

I if \u\ > 2n. 

We consider then the family (^") neN of potentials over R given by 

1 2 

§ n ( u ) = o (Au,u) - T G n {u d ). 
1 n 

This family satisfies the following properties: 

(i) For every n S N the potential (j)" is of class C 2 and b" = — V(|)" is globally Lipschitz. 

(ii) For n < m G N we have b n = b'" over the region IT' = {u G R d : \ud\ <n}. 

A/1 

(iii) For every n G N we have lim inf|„|^. +00 > 0. 

Since b" is globally Lipschitz, for each u£R d there exists a unique solution to the ordinary differential 
equation 

U n ' u = b n (U n ' u ) 
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with initial condition u. Such solution is globally defined and describes the same trajectory as the 
solution to (11.11 ) starting at u until the escape from II". In the same way, for each x £ W l and £ > 
there exists a unique global solution to the stochastic differential equation 

(5.1) dU"' u ' e = b"(U n '"' e )dt + edW 
with initial condition u. 

As before we use [/"'" for U n ' u '°. Since b" coincides with b over the ball B n (0) of radius n centered 
at the origin, if we write 

x?" = mf{t > : \U n ' u ' e (t)\ > n}, x" := lim x£' M , 

then for ? < Xg we have that U u ^{t) := lim„^ +00 [/"'"' e (f) is a solution to 

(5.2) dU u ' e = b(U"' e )dt + edW 

until the explosion time x" with initial condition m. Moreover, if we define the stopping times 

j#" =inf {?>(): £/"'"• E (0 ^IT 1 }, 
it can be seen that (ii) implies that 

x e = hm 7i e 

n— >+<*> 

and that U u ' e coincides with the process U n,u ' e until the escape from II". On the other hand, (i) 
guarantees that for each n G N and u £~M. d the family (^"'"' E ) E>0 satisfies a large deviations principle. 
Finally, from (iii) we get that there is an unique invariant probability measure for the process U n,s for 
each £ > given by the formula 

^£ J A 

WM 

where Zg = J R de e 2 du. Hereafter, when we refer to the solution of (15.21) we mean the solution 

constructed in this particular way. 

5.2. Proof of Theorem 12.41 We split the proof of Theorem l2.4l in two parts, the first one is immediate 
from the continuity of the solutions of (11.21 ) with respect to £ in intervals where the deterministic 
solution is bounded. 

Proposition 5.1. For any fixed 8 > and u € D e we have 

limP„(Xg <x -8) =0. 

e->0 

Proof. We may assume that x(j > 8 since the proof is trivial otherwise. Now, as the deterministic 
system U" is defined up until Xq, if we take M := sup 0< , <T «_ g \U" \ < +°° then x" < X M — 8 implies that 

sup \U 2M ' l '- £ {t)-U 2MM (t)\ > 1. 
0<t<%,-8 

By (|II]> we get (EB- □ 
Proposition 5.2. For any 8 > and u € D e we have 

limP„(x e > x + 8) =0. 

Moreover, the convergence is uniform over compact subsets of D e . 
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Proof. Fix 8 > 0, X a compact set contained in D e and let Y" be the solution to the ordinary differential 
equation 

2 



yll 



AY U + -Y»' b e d 



with initial condition u € By the Comparison principle we have that U" > Y" for as long as U" is 
defined. Since Y" is the solution to a linear system of ordinary differential equations whose associated 
matrix is symmetric and negative definite, we get that there exists G R such that for all u G %, 
every coordinate of U" remains bounded from below by + 1 up until Xg. If for p G R and M > 
we write 

0^ := {y G R d : y k > p for all < k < d - 1 , y d > M} 

then by the Maximum Principle and the previous statement we have that T u := inf {t > : U" G } 
is finite. Moreover, as U M+ agrees with U" until the escape from Hm+2, we obtain the expression 
T u = mf{t > : tjM +2 ' u £ 0^+* }. Taking T x := sup ue5C T u < +oo we may compute 

p„(x e (0fj > t u ) <p I( (7if +2 Ax e (0fj > t u ) +p l( (7if +2 < r M ,x e (<) > r B ) 



<2P„ sup \U M+2 > e (t)-U M+2 {t)\ > 1 
^0<f<7; 



<2P !( sup |t/ 



M+2,E/ 



■U M+2 (t)\ > 1 



0<?<T y 



from which by (12.21 ) we obtain 
(5.3) 



limsu P P ii (x e (0fj >T U ) =0. 



uexi 



On the other hand, by the strong Markov property for U u ' e we get 

P u (x E > xo + 5) < P„ (x E > P„ + 5) < sup P v (x E > 5) + sup P u (x e (0f ) > P„) . 



Taking into consideration (15.3I ). in order to finish the proof we only need to show that the first term on 
the right hand side tends to zero as £ — > for an adequate choice of M. To see this we consider for 
each £ > and y G &^ the processes Y y,e and Z V E defined by 



dY 



y,E 



2 



AY y ' E + -Fj' b e d ) * + £^W, 



and Z v £ := C/ y ' e — Y y,e , respectively. Notice that since Y y,e is globally defined and both U y ' e and Z y 
have the same explosion time. Also note that Z v £ satisfies the random differential equation 



dZ 



y,e 



AZ y >* + - 



V,E\ + 



7 V,E 



The continuity of trajectories allows us to use the Fundamental Theorem of Calculus to show that 
almost surely Z V E ((o) is a solution to the ordinary differential equation 



(5.4) 



(0(o>) = -AZ» e (co) + \ ( [(t/f ) + ] "((0) -Zf (CD))* 
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For each y G and £ > let £>g be a set of probability one in which ( 15.41 ) holds. Notice that for 
every co G Qg we have the inequality 

Z-^)>-AZ>''>)-^Zf (<Q)e d . 

Using the Comparison Principle we conclude that Z- v,E (co) > for every CO G £lg and, therefore, that 
the inequality U y ' e ((0) > Y ye ((i)) holds for as long as U y ' e ((i)) is defined. 

For each y G &^ and £ > let us also consider the set 

Q y e = \(0£Q : sup |y- v ' E ((0,?)- Y y ((i),t)\ < 1, sup |eW(a>,f)| < l}- 

0<f<8 0<f<8 J 

Note that lim e ^o inf ye @M P(Cl y e ) = 1. Our goal is to show that if M is chosen adequately then for fixed 
y G ©j!f the trajectory £P' e ((q) explodes before time 8 for all co G £> e PI Hg. From this we get that 

inf P(&1)= inf P(Q&n££) < inf P v (x e < 8). 

and by letting E^Owe conclude the result. 

So let us take y € ©jf , CO £ fig PI f2 e and suppose that U y ' e ((0) is defined in the interval [0, 8]. Notice 
that since CO E ilg D & e then the (d — l)-th coordinate of F V E (co,f) is bounded from below by — 1 
for t € [0,8]. By comparison we know that the (d — l)-th coordinate of £//' E (co,f) is bounded from 
below by — 1 as well. 

From here we deduce that the last coordinate of £P' E (co) verifies the integral equation 

Uf(a,t)>U>;\(»,s)+j'^(-Uf(«),r)+p x -l+hg(U y / 

for s < t G [0, 8] . We can take M G N large enough to guarantee that there exists a constant a > such 
that for all m>Mwe have 

2 

~2 ( - m + p x - 1 + hg(m)) > am p . 
If we recall that £/J' E (co,0) > M then our selection of M implies that 

U-;\(o,t) >M-\+aj^ (Uj\(o,u)) p du 

for all t G [0,8]. But if this inequality holds and M is large enough, one can check that [/^(co) 
explodes before time 8, which contradicts our assumptions. Therefore, if y G ©j!f and co G £l e D ^ E 
then £/- v ' E (co) explodes before time 8 and this fact concludes our proof. □ 

Combining these two propositions we get Theorem |2.4j Observe that the bounds obtained decay to 
zero exponentially fast due to Proposition ( I2.2I ). 

6. METASTABLE BEHAVIOR FOR INITIAL CONDITIONS IN D 

Finally we focus on initial data in Do, where the metastability phenomenon can be appreciated. We 
start with the construction of an auxiliary domain that contains the origin and such that the exit time 
from this domain is asymptotically equivalent to the explosion time. 
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6.1. Construction of an auxiliary domain. In order to proceed with our analysis of the explosion 
time we must first construct an auxiliary bounded domain. The puipose behind this construction is 
to reduce our problem to a simpler one, the escape from this domain. This is easier because we may 
assume that the drift coefficient b is globally Lipschitz, as the escape only depends on the behavior 
of the system while it remains inside a bounded region. In this case, large deviations estimates as 
the ones proved by Freidlin and Wentzell apply. We need a bounded domain G which verifies the 
following properties: 

(1) G is bounded, contains 1 and the origin. 

(2) There exists c > such that fi c (0) C G and for all y G B c (0) the system U y is globally defined 
and tends to zero without escaping G. 

(3) The border of G can be decomposed in two parts: 3 1 and dG \d l . The region of the border 3 1 
is closed and satisfies mm Me 9 G (|)(w) = min 1(£ gi §(u) and 

inf fo(u) > min(b(j<). 

uedG\d' uedG 

(4) For all y G d l the deterministic system U y explodes in finite time. 

The domain G can be constructed as follows. Let us consider the value of (j) at the saddle point 1, 
(j)(l) = -1/(^ + 1) + 1/2 >0 = (j)(0) and c> such that (j)(w) < (j)(l) for u G5 C (0). 

For each point u G 3fi c (0) consider the ray R u := {ku : A > 0}. Since the vector 1 is not tangent 
to %>[ at 1, we may take a sufficiently small neighborhood V of cl such that for all u G V n3fi c (0) 
the ray R u intersects W{ D (M>o) d - For such V we may then define A„ = inf{X > : Xu G W/} for 
ueVf] dB c (0). If we consideiQ 

T) := inf §(l u u) > <|>(1) 

»e3[yn3B c (0)] 

then the fact that (|)(t/ (t)) is strictly decreasing (see Proposition 13 .11 ) allows us to shrink V into a smaller 
neighborhood V* of cl such that (j)(v) = r| for all v G n3fi c (0)]. Let us also observe that since 1 is 
the only saddle point we can take V sufficiently small so as to guarantee that max{§(ku) : X > 0} > r\ 
for all u £ dB c (0) \ V*. Then if we take the level curve = {x G W 1 : §(x) = r\} every ray R u with 
u G dB c (0) \ V* intersects C^. With this we may define for each u G dB c (0) 



X l( if u G V* 

inf{X>0:XMGC T1 } if u G B c (0) \V* 



Notice that the application u i->- A,* is continuous. Due to this fact, if (5 := {Aw : < A < X* , m G dB c (0)} 
then 3(5 = : u G 3fi c (0)}. To finish the construction of our domain we must make a slight radial 
expansion of G, i.e., for a > consider G defined by the formula 

G:={Xu:0<X< (1 + a)l* u , u G dB c (0)}. 

Let us observe that Theorem 13.51 insures that G verifies condition (1). Since A,* > 1 for all u G dB c (0) 
then it must also verify (2). Also, if we define 3 1 := {(1 +a)A*(w) : X*(u)u G V*} then 3 1 is closed 
and if a > is taken small enough then (3) holds. Finally, due to Theorem 13 .5 1 we have 3 1 C D e and 
so (4) is verified. See Figure I67TI 



By d[V ndB c (0)] we mean the border of the (d — 1) -dimensional manifold V n3B c (0). 
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Figure 2. The level curve CV, and the stable manifold of 1. 



6.2. The escape from G. The behavior of the explosion time for initial data u € Do is proved by 
showing that, with overwhelming probability as £ — s> 0, the stochastic system describes the following 
path: 

(i) It enters a neighborhood of the origin B c (0) in before a finite time T that does not depend on 
£. 

(ii) Once in B c (0) the system remains in G for a time of order e A / e and then escapes from G 
through 3 1 since the barrier imposed by the potential is the lowest there. 

(iii) After escaping G through 3 1 the system explodes before a finite time x which does not depend 
on £. 

The fact that the domain G is bounded allows us to assume that b is globally Lipschitz if we wish 
to study the behavior of our system while it remains inside G. Indeed, we may take jiq € N such that 
G C B„ (0) and study the behavior of the solution to (15.11 ) since it coincides with our process until 
the escape from G. Then we can proceed as in the double-well potential case to obtain the following 
results (see iflOl pp 295-300] for their proofs). Hereafter, B c (0) denotes the neighborhood of the 
origin highlighted in the construction of G in the previous section. 

Theorem 6.1. Given 8 > we have 




Theorem 6.2. The stochastic system verifies 



lim sup P u (U e (? e (dG)) <£d l ) = 0. 



ueB c (0) 



From these two theorems we can obtain the following useful corollary. 



Corollary 6.3. For any 8 > we have 



lim sup P u [% E (d )>e~S r )=0. 
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□ 



Proof. One can easily check that 

sup P u (x e {d l ) >e^) < sup P u U e (dG)>e^) + sup ^fe e(9G) g 9 1 

u£B c (0) V ' ueB c (Q) V 7 h£B c (0) 

Concerning the asymptotic distribution of x e (3G) we can obtain the following result. 
Theorem 6.4. Let y e > be defined by the relation 

P (x & (dG)>ye)=e- 1 . 
Then there exists p > such that for all t > we have 

lim sup |P„(x e (3G) > /Ye) - = 0. 

e ^°«GBp(0) 

6.3. Bounds for the explosion time. This section is devoted to the lower and upper bounds for the 
explosion time. More precisely, in this section we show that given 8 > 0, for all u G D one has 

limP H (x e < e^) = 

£^0 V / 

and 

/ A+5\ 

hmP s T £ >e7 =0, 

E^O V / 

where the convergence can be taken uniform over compact subsets of Do- The proofs of these bounds 
essentially follow ifTOl . where analogous bounds are given for the tunneling time. However, unlike 
the double-well potential model, the use of localization techniques becomes necessary at some points 
throughout our work. We begin first with the lower bound. 

Proposition 6.5. Given 8 > and u G Do we have 

(6.1) \imP u U e <e^) =0. 

Moreover, the convergence is uniform over compact subsets of Do. 

Proof. First observe that since for u G G we have P u (x E > x e (3G)) = 1 then (16.11 ) holds uniformly over 
any small neighborhood of the origin by Lemma loTTI Next, we generalize the result for any u G Do- 
For each u G Do there exist T u > 0, 8„ > and n u G N such that the deterministic system beginning at 
u reaches 5j> (0) before T u , remaining in B Hu (0) and at a distance 8„ from dB„ u (0) on [0, T u \. It follows 
that £/"'"" does so as well. From this we obtain 

P„(T e (5p(0)) > T u ) < P u (mm{z'i" , X E (fi p (0))} > T u ) +P U {^" < T u ) 



<P U ( sup \U"- £ (t)-U n "(t)\>^)+pJ sup \U n ^(t)-U n "(t)\>^). 
y o<t<T„ 2 / ^o<«7; 2 / 



Using estimation (12.21 ) for the family (^""'"' E ) e> o we conclude 
(6.2) UmP u (T e (fl p (0)) >T U ) =0. 

Therefore, if we write 

P u [% < e^*) < P u (x e (B p (0)) < x E < e^) +P M (x e < T u ) + P l( (x E (fi p (0)) > T u ), 
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then the last two terms on the right tend to zero when £ — > as a consequence of what we stated above. 
By the strong Markov property for £/" ,e we have 

P„(x e (fl p (0)) <x £ <e^) < sup P y (%t<e*r)< sup P v (x"°(dG) < e^) 

V ' V6B P (0) V ' yeBp(O) V ' 

where no is taken as in the first step. Since the rightmost term tends to zero by Lemma loTTl we conclude 
the result for arbitrary u € Dq. The uniform convergence over compact subsets %, of D is proved in 
a similar fashion by taking S i( and T u uniformly over OC as in Proposition 15.21 □ 

Now we turn to the proof of the upper bound. As we stated before, when studying the behavior of 
the stochastic system under initial conditions u € G and for small £ > we typically observe that the 
process U" ,e escapes from G through 3 1 since the cost imposed by the potential is the lowest there. 
Once in 3 1 the influence of noise becomes negligible and the process then describes a path similar 
to the deterministic trajectory until exploding in a finite time. We formalize this statement in the 
following proposition. 

Proposition 6.6. There exists To > such that 

limsupP„(T E > 7b) =0. 



£^0 



ue3 



Proof. Since 3 1 is a compact set contained in D e , the proof follows from Proposition 15.21 and the fact 
that sup„ e91 x° < +oo □ 

With this proposition we are able to conclude the upper bound. 

Proposition 6.7. For each 8 > and u G Do we have 

/ A + 8 \ 
limP i( x e > = 0. 

e^o V / 

Moreover, the convergence is uniform over compact subsets of Do. 

Proof. We proceed in two steps. 
1. We check that given 8 > we get 

li 



(6.3) lim sup P x x E > e~&~ = 0. 

It is not hard to show that for £ > small enough the strong Markov property yields 

sup P u (x e >e^] < sup pJ x E (3 1 ) >e^ J +supP i( (x e > 7b) + sup P«(t/£( 9G) £ 3 1 
wefi c (0) »eB c (0) \ / ued 1 ueB c (o) 



where To > is taken as in Proposition 16.61 We finish this first step by observing that the right hand 
side converges to zero. Indeed, the first term does so by Corollary 16.31 the second by Proposition 
and the third by Lemma [6T2l 

2. We now generalize the result for u € Dq. This follows from the fact that 

Puf'h > e^) < sup PJ x e > tL. ) +P„(x e (B c (0)) > T u ) 
v ' ueB c (0) V 1 ) 
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by the strong Markov property. Observing that the first term on the right hand side of the equation 
tends to zero by (16.31 ) and that the second term does by (16.21 ). we obtain our result. The convergence 
over compact subsets of D can be seen as in Proposition 15.21 □ 

6.4. Asymptotic distribution of the explosion time. Our main objective in this section is to prove 
the asymptotic memory loss of the normalized explosion time ^ . The proof focuses on studying the 
escape from G. The asymptotic memory loss for x E can be deduced once we show that the time in 
which the process exits from G and the explosion time are asymptotically similar. We formalize this 
last statement in the following proposition. 

Proposition 6.8. There exists a positive constant To such that for all u € Do fl G 

limP t (x E >x e (3G) + ro)=0. 

e->0 

Proof. Let us observe that by the strong Markov property 

P„(x e > x e (9G) + T ) < supP v (x e > T ) + P„(x E (3G) < x e (B c (0))) + sup P ( (£/ T e e(9G) £ d l ) . 

yed 1 ueB c (0) 

We can now conclude our desired result by the use of Proposition 16. 61 and Lemma I6T21 □ 

We are now ready to establish the asymptotic memory loss of the explosion time. Having the former 
proposition at our disposal, the rest of the proof is very similar to the one offered in the double-well 
potential model. We emphasize that the main difference with this case lies in how to show this last 
proposition. In the double- well potential the corresponding statement to Proposition 16.81 holds due 
to the fact that the tunneling time for initial conditions in the deepest well is of order one. This can 
be easily deduced from the Freidlin-Wentzell estimates. Analogously, in our model Proposition 16.81 
holds since now the explosion time for initial data in D e is of order one. However, the lack of a global 
Lipschitz condition forces us to proceed differently in order to show this last fact. We recall that a 
proof of this is contained essentially in Proposition 15.21 We now give a brief sketch of the rest of the 
proof of Theorem 12. 61 in the following lines and refer to for further details. 

Sketch of proof of Theorem 1276] 

(1) We first check that, for p > small enough, lim e ^o su P l , e 5 p (o) \Pu{te(dG) > tfi e ) — e~'\ = 0. 
This is due to the fact that lim e _ ! .o y = 1 . 

(2) Next, we prove that Po(x e > ?Pe) = e ~' for t > 0. This is done with the help of Proposition 16. 8 1 
and the previous step. 

(3) With the help of appropriate coupling techniques we establish the uniform convergence over 
any small enough neighborhood of the origin. 

(4) Finally, by using the strong Markov property, we conclude the result for arbitrary initial data 
u € Do- 
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